Charge and spin collective modes in a quasi-lD model of Sr 2 Ru04 



Suk Bum Chung 1 , S. Raghu 1,2 , Aharon Kapitulnik 1 , and Steven A. Kivclson 1 

1 Department of Physics, Stanford University, Stanford, California 94305, USA and 
2 Stanford Institute for Materials and Energy Science, 
SLAG National Accelerator Laboratory, Menlo Park, California, 94025, USA 

(Dated: March 8, 2013) 

Given that S^RuCU is a two-component p-wave superconductor, there exists the possibility of 
well defined collective modes corresponding to fluctuations of the relative phase and spin-orientation 
of the two components of the order parameter. We demonstrate that at temperatures much below 
T c , these modes have energies small compared to the pairing gap scale if the superconductivity 
arises primarily from the quasi ID {d xz and d yz ) bands, while it is known that their energies become 
comparable to the pairing gap scale if there is a substantial involvement of the quasi 2D (d xy ) band. 
Therefore, the orbital origin of the superconductivity can be determined by measuring the energies 
of these collective modes. 



I. INTRODUCTION 

The layered pcrovskitc material Sr2Ru04 has attracted 
strong interest for over a decade because of experiments 
indicating that it is a spin-triplet superconductor— ~— with 
spontaneously broken time-reversal symmetry£>£. The 
only state that possesses both of these properties in a 
tetragonal system with spin-orbit coupling is the chiral 
p-wave state, the electronic analogue of 3 He-Ai£. In 
a single-band, quasi-two dimensional system, this state 
is expected to be a topological superconductor: it is 
fully gapped and has topologically protected Majorana 
fermion zero modes in vortex cores and along its edge. 

However, the observation of power laws in specific 
heali 9 - and NMRji£, the absence of electric currents along 
edges and domain wallsii, and the absence of a split 
transition in the presence of an in-plane magnetic fielcU^ 
are sharply inconsistent with the theoretically expected 
properties of a simple chiral superconductor. Motivated 
by these inconsistencies, three of usi 3 . proposed that the 
multi-band nature of the material is essential. 

The electronic structure of this system is derived from 
the Ru t2 g electrons {d xz ,d yz , and d xy ). These orbitals 
produce three Fermi surfaces, denoted a, /3, and 7. The 
a and ft surfaces are derived from the {d xz , d yz } orbitals 
and are quasi-one dimensional (ID), whereas the 7 sur- 
face is derived primarily from the d xy orbital and is quasi- 
two dimensional (2D) (see Fig. []}. Because of the strong 
differences in character between the different bands, su- 
perconductivity in this system is likely derived primarily 
from either the d xy , or {d xz ,d yz } orbitals: in either sce- 
nario the "active" electrons induce superconductivity in 
the remaining "passive" subset via the proximity effecli^. 
However, to the extent that the proximity effect is weak, 
there is a range of circumstances in which the exper- 
imental observations reflect superconducting properties 
of mainly the active orbitals. 

There are sharp distinctions between the two possibil- 
ities for the active electrons. For example, the p x + ip y 
state obtained only from the a, j3 bands is topologically 
trivial because the a and ft bands form, respectively, a 
hole and electron pocket leading to a net zero Chern 



number 13 . If these were the active electrons, the sys- 
tem behaves as a topologically trivial system near T c . 
At T = 0, superconductivity occurs on all three bands 
and is therefore topologically non-trivial; however, the 
size of the induced gap on the passive 7 Fermi surface 
can be substantially smaller than those on the a, /3 sur- 
faces, making it difficult to experimentally detect the 
topologically non-trivial character of the ground state. 
This scenario could explain, for instance, the absence of 
any detectable edge currents in the systemii. By con- 
trast, when 7 is the active band, the system behaves as 
a topological superconductor even near T c . Recently, it 
was shown that there is an intrinsic Kerr response near 
T c only when {a, 0} are the active band a 15 i 16 . In this 
paper, as part of a further exploration of the experimen- 
tally accessible properties that could distinguish between 
the two cases, we study the qualitative differences in the 
character of the charge and spin collective modes. 

In forming the a and /3 bands, the d xz and d yz or- 
bitals are coupled to each other only via spin-orbit cou- 
pling and second-neighbor hopping terms, both of which 
are relatively weakii in S^RuO,}. The limit in which 
these orbital mixing terms vanish defines a multicritical 
point at which the superconducting state breaks a higher 
[5'0(3) S pi n x C/(l)chargo] 2 symmetry. Proximate to this 
multicritical point, i.e. for weak orbital mixing, there are 
low energy "almost Goldstone modes" associated with 
fluctuations of the relative phase and spin orientation 
of the x and y components of the superconducting or- 
der parameter. By contrast, the same collective modes 
hav o 18 ' 19 energies of order Ao when 7 is the active band. 

This paper is organized as follows. In section II, we 
discuss the general physics of the collective modes in a 
multi-component superconductor at a qualitative level 
and describe the general form of the non-linear sigma 
model (NLSM) of the system. In sections III-V, we de- 
rive the NLSM from the microscopic physics of the quasi- 
lD model and obtain the gaps of the collective modes. 
Lastly, we discuss schemes for detecting the modes and 
consider their broader implication for the superconduct- 
ing properties of the system. 




FIG. 1. Schematic view of the Ru d xz and d yz orbitals in 
the Ru-0 plane. The wave-function overlap is larger along 
the direction of the black arrows than along the direction of 
the gray arrows, giving us t > t' , and therefore, the quasi-lD 
nature of the bands originating from these orbitals. 



II. COLLECTIVE MODES IN 
MULTI-COMPONENT SUPERCONDUCTORS 



A. Qualitative discussion 



A superconductor described by a multi-component or- 
der parameter (p x and p y in the present context), will 
have collective mode excitations associated with the rel- 
ative phase difference, <f>- = 9 X — 9 y . At zero tempera- 
ture, such a mode would be expected to have a frequency 
fj^o ~ \JJ lx, where \ ~ -W(O) is the compressibility 
(A(0) is the density of states at the Fermi energy) and 
J is the second derivative of the condensation energy 
with respect to (/>_. Given that the condensation en- 
ergy ~ iV(0)|Ao| 2 where |Arj| is the root mean squared 
gap magnitude, this generally means that fiwrj ~ l^o|- 
Similar considerations apply to fluctuations in the rela- 
tive orientation of the spins (i.e. the d- vector) in a two- 
component triplet superconductor. 

By contrast, if the two components of the order pa- 
rameter are associated with different orbitals, i.e. the 
p x component with the d xz and the p y with d yz or- 
bital respectively, and if mixing between the different 
orbitals were absent, then collective fluctuations of the 
relative phase and spin-orientation would be gapless. The 
orbital-mixing terms, which we schematically denote as 
SH (and which we discuss in detail below), result in a 
non- vanishing dependence of the condensation energy on 
4>-. As a result, if the superconductivity arises primarily 
from the quasi-lD bands, then the relative phase mode is 
expected to have an energy tuvo ~ 7A0, where 7 vanishes 
continuously as 8H tends to zero. Naturally, similar con- 
siderations apply to the relative spin orientational fluc- 
tuations. As we explain below, since appreciable orbital 
mixing occurs only where the bands cross one another, 
7< 1 even when the characteristic scale associated with 
8H^> |A |. 



B. Non-linear sigma model 

To analyze the low-lying collective modes in a two- 
band spin-triplet superconductor we consider the non- 
linear sigma model (NLSM) valid deep inside the super- 
conducting phase. We express the order parameter as 

A a;ss ,(fc) = Af a QL)e ie "(ia 2 d a ■ a) ss ,, (1) 

where a = x, y labels the two components of the or- 
der parameter, s, s' the spin indices of the two electrons 
forming the Cooper pair, A the order parameter ampli- 
tude, / Q (k) the pair wave-function which is determined 
by the microscopic form of the pairing interaction, and 
d Q a real unit vector in spin-space. In an A-orbital ba- 
sis, f a (k) (and consequently A Qj g a /) are N x N matrices 
which transform as the x and y components of a vector 
of the tetragonal point group. 

The NLSM action is obtained by holding A and f a (k) 
fixed, and focusing on the long-wave-length fluctuations 
of 9 a and d a : 

C= \Y. [ w I W 2 + M(dtd a ) ■ (d t d a ) + ... 

a 

—J[6x — @y, dx ■ d y ] + . . . 

-^r (d a -d a /3 -«) + ••• (2) 

a 

Here, the first line contains the terms which respect the 
symmetry 

of the multicritical system. The " . . . " includes terms 
proportional to spatial-derivatives and higher powers of 
time derivatives. Our interest is in the long wavelength 
limit and we shall neglect spatial derivatves; however, 
these terms arc the only ones which reflect the underlying 
tetragonal lattice symmetry. 

The second-line includes terms derived from inter- 
orbital mixing without spin-orbit coupling, and the third 
line includes spin-orbit coupling terms. Again, deriva- 
tive terms are neglected. Here, J[(f>,x\ = J[<p + ir, x] = 
J\— <f>, x] = J[<j>,—x]. Note that, there is no Goldstone 
mode resulting from fluctuations of the overall phase 
<f> + = 6 XZ + 9 yz due to the Andcrson-Higgs mechanism. 

In the next section, we define a microscopic model of 
the superconductivity on the quasi-lD bands and from 
it derive estimates of the various couplings that appear 
in Eq. [2] Estimates for W and AA, both of which are of 
order l/iV(0), are discussed below Eq. [5J To first order 
in the (spin rotationally invariant) orbital mixing terms, 
the form J is found in Eqs. [T51 and [TBI to be 

J[4>-J ■ d'} - Jo cos(20_) [2(d ■ d') 2 - l] (3) 

and an estimate of Jq is presented in Eq. 1171 Finally, 
an expression for Tq to second order in the spin-orbit 
coupling is given in Eq. 1191 

Because J7o > 0, the orbital mixing term is minimized 
either when d x is parallel to d y and </>_ = 7r/2, which cor- 
responds to the chiral superconducting state analogous to 
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the A- phase of 3 He, or when d x is perpendicular to d y 
and (f> = or 7T, which corresponds to the time-reversal in- 
variant superconducting analogue of the B-phasc of 3 He. 
The degeneracy between these two phases is exact (at 
mean-field level) in the absence of spin-orbit coupling; 
this observation and the various periodicity conditions 
impose general constraints on the form of J , which are 
satisfied by our present result. The degeneracy between 
the A and B phases is lifted by spin-orbit coupling: for 
Tq > 0(To < 0), the A(B) phases have larger condensa- 
tion energy. 

Wc now proceed to show explicitly how these quali- 
tative considerations apply in a simple, but physically 
motivated microscopic model of Sr2Ru04. 

III. QUASI-1D MODEL OF SR 2 RU0 4 
SUPERCONDUCTIVITY 

As a microscopic representation of the problem of the 
pairing in the quasi-lD bands, we consider an idealized 
form of the Bougoliubov-de Gennes Hamiltonian for the 
quasi-particles in a p-wave superconducting state, 

HsdG = H mu ui + SH. (4) 

Here H mu ui represents the multicritical point model, in 
which there is no orbital mixing or spin-orbit coupling, 
and SH (assumed small) represents terms which break 
the higher symmetry of the multicritcal point. 

A. The multicritical point Hamiltonian 

Hmuiu is the mean-field Bogoliubov-de Gennes Hamil- 
tonian for decoupled xz and yz orbitals: 

H mu lu= ^2 X^ a ( k ) C aks C aks (5) 

a—xz,yz ks 

+\ E E E[ A ^( fc )i k X- k v + h - c -]' 

a—xz,yz a—x.y k;ss' 

where c a k s is the annihilation operator for an electron 
with the momentum k, orbital index a — xz,yz and 
spin polarization s, £ a (k) = e a k — [A, the chemical poten- 
tial /i is set so that these bands are two-third fille d 17 ' 20 , 
and A^ a ^, is the appropriate orbital diagonal matrix 
element of the pair-field defined in Eq. [TJ We fur- 
ther simplify the model by taking the band-structure 
in the absence of orbital mixing to be strictly one di- 
mensional, e a k = — 2tcosk a and the x and y com- 
ponents of the order parameter to originate entirely 
on the corresponding ID band, A^f], = A^' g , = 0, 

Ag,(fc) = Asm(k x )e e *[ia 2 d x ■ a] ss , and A^,(fc) = 

Asm.(ky)e ey [i<j2d y ■ a] ss '. While this simplified band 
and gap structure simplifies the explicit calculations, the 
qualitative results we have obtained are not affected by 



the inclusion of moderate transverse (but still orbital- 
diagonal) components of the hopping matrix. 

Because the two bands are related by rotation by n/2, 
H mu iti respects the tetragonal symmetry of the material 
even though each band is one dimensional. From the 
fact that A x lives entirely on the d xz band and A y on 
the d yz band, it further follows that H mu i t i respects a 
full [5 f O(3) sp i n x C/(l)chargc] 2 symmetry, so there is no 
dependence of the free energy on 9 a or d a . In particular, 
note that the quasiparticle spectrum is fully gapped not 
only in the A-phase (e.g. for d x -d y = 1 and 6 x —6 y = it/2) 
and B-phase (e.g. for d x • d y = and 6 X — 9 y = 0), but 
also for a p x ± y (e.g. d x ■ d y = 1 and 8 X — 8 y = 0, tt) which 
might otherwise have been expected to posses gap nodes. 

B. Orbital mixing and spin-orbit coupling 

SH contains all permissible terms in the quasi-particlc 
Hamiltonian which break the [S , 0(3) sp i n x C/(l) c harge] 2 
symmetry of the multicritical point, of which the most 
important are band-structure terms that mix the two or- 
bitals, 

^ = E Ak ( c kk, s c ^,k, s + h.c.) (6) 

+ * E E e °b a ss' c Ls c b^ ,+■■■ 

a : b—xz,yz k;ss' 

where the first term reflects second-neighbor hopping the 
xz and yz orbitals, so Ak = 25tsa\k x smk y , and the sec- 
ond term represents the Ru atomic spin-orbit coupling 
where i c ab = ie a bc are the spin-1 matrices representing 
the effective orbital angular momentum of the tig or- 
bitals and er a are the usual Pauli spin-matrices. We can 
see that if we include the xy orbital, the atomic spin- 
orbit term will include the additional terms £ x a x +£ v <j v . 
These terms, however, do not alter any of our qualita- 
tive results, except those we discuss in Section ITVCl and 
Appendix [B] 

A more detailed model of the electronic structure of 
Si'2Ru04 might include additional terms in SH (indicated 
by . . .) but we will show that these terms are certainly 
smaller than the terms we have kept and therefore do 
not qualitatively effect the outcome of the calculations, 
so long as we focus on the limit 

t»l<yt|~H»|A|. (?) 

Eq.© omitted the intcrobital interaction terms of the 
form 

SH int = V d iz,,a d lz^' d yz^' d *z,ia ( 8 ) 
iaa' 

that are local and therefore play a role everywhere in mo- 
mentum space. However, in the BCS ground state, such 
interactions cause scattering only in the Cooper channel, 
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and therefore play an appreciable role only in momen- 
tum points where the bands cross. Additional interac- 
tion terms such as the interorbital singlet-pair hopping 
terms do not play a role in the spin-triplet supercon- 
ductor. Among higher order interaction processes, only 
those terms proportional to the interorbital susceptibility 
(vanishingly small in the present case) affect the interor- 
bital modes. Thus, we are led to the same conclusion: 
the mixing between the different orbitals whether it be 
due to tunneling or to interactions, is weak. 



IV. DERIVATION OF THE NLSM 

In this section, we will compute the terms in the NLSM 
in Eq(5]from the microscopic model in Eqs. 0K1 

A. Kinetic terms, W and M 

Firstly, we derive expressions for the time-derivative 
terms in Eqj2] which do not involve mixing of the orbitals 
or spin-orbit coupling; to leading order, these terms can 
be computed in the limit SH = and A = 0. Because 
the number density is canonically conjugate to the phase 
of the superconducting order, these terms ultimately re- 
flect the energy cost of shifting electrons from one orbital 
to another or from one spin polarizing to the other. To 
quadratic order, the cost in energy per unit volume asso- 
ciated with a transfer of charge Sn per unit cell from one 
band to the other and with a spin density S a in band a 
is 



K[6n,S a ] = -^-(Sn) 2 + -^- 

4Xc 2 Xs 



(9) 



where Xc = N(0) and Xsp — AT(0)/4 are the density 
and spin susceptibilities for each quasi-lD band. In the 
standard fashion, spin-triplet Cooper pairing gives rises 
to the commutation relation e! 2 1-23 



[Sn, e H 



-Me 



i<j>- 



[sidi]=i^ k 5 ab dt, (io) 

thus, we can regard K as the 'kinetic' energy density of 
4>— and d a . By using the Heiscnbcrg equations of motions 
dt4> = i>[4 > >K]/fo an d dtd a = i[d a , K]/H, we obtain 



K[<t>X] = 



16 



Xsp 



E(w 2 . (ii) 



This is equivalent to the first line of the NLSM action of 
EqG3 



B. Relative Josephson coupling, J 

The leading contribution to the interorbital Josephson 
coupling, J , can be computed in the n = limit. As a 



warm up, we first calculate the ground state energy for 
the spinless case where the relative phase between the 
condensate of two components of the order parameter is 
set to (f>_ = 6 X — By. The ground state energy can be 
computed as 



£bfa-) = -53^«*(0-)-eak]/2, 



(12) 



where £^ a k(<^-)'s are the eigenenergies of the spin "up" 
quasiparticles in T-L BdG i n Eq. [5] with d x = d y = z and 
A = 0. As we show in Appendix I A 11 the local stability of 
the chiral state follows from from the fact that £'q{4>- = 
±tt/2) = and £' Q ' (0_ = ±tt/2) > 0. Indeed, since in 
the limit of weak mixing, we know the dependence of 
the energy on 0_ must be approximately harmonic, it 
follows that Eq is related to the value of the Josephson 
coupling. Time reversal symmetry implies that £o (</>-) = 
£o(—<f>-) and the fact that A a transforms like a vector 
under rotations by n/2 implies that £o (</>—) = £q(tt— 4>-). 
Thus, 



£o(4>-) ~ const. + Jo cos(20_), 



where 



Jo - -jEq 



ok 



(13) 



(14) 



We have computed Jq by numerically diagonalizing the 
model in Eq. |4] with r) = - see Fig. [2] and Appendix 

ixn 

It is straightforward to generalize the above result to 
the spinful case and obtain J of Eq.Q. From the defi- 
nition of the d- vecto r 20 ! 23 



It 



A/ Q (k)e i9 °d 



(15) 

it follows that it is always possible to chose the spin 
quantization axis, e z , to be perpendicular to both d x 
and d y , in which case all pairing is between like spins 
(Aq-j-j, = 0). Thus, for r\ = 0, the Josephson coupling is a 
sum of terms two equal contributions from spin-up and 
spin-down pairs: 



£ o (0tt ) = Jo (cos 2(j)^ + cos 2<%) 



where 



Hi 



(16) 



and 



?tt — °x - Vy -r a -i- 
cos a — d x ■ dy. A small exercise in trigonometry thus 
leads to from this result to Eq. |3]with Eq. [14] for Jq. 

We have numerically evaluated the integral over k in 
Eq. [TUto obtain the results for Jq shown in FigJ5J It is 
clear that Jq scales as 



Jo 



\St\ A 2 \St\ 



r 2 



t 



A(0)A 2 , 



(17) 



where N(0) — Wir/ (v3i) is the density of states for each 
quasi-lD orbital; this can be also obtain through analytic 
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FIG. 2. The left panel shows as the black solid curves a 
schematic plot of the Fermi surface originating from the quasi- 
1D bands which derive primarily from the Ru d xz and d yz 
orbitals; it was computed from the the microscopic model of 
Section II. (The blue dotted curve represents the Fermi sur- 
face of the quasi 2D band arising from the Ru d xy orbitals, 
which are not included explicitly in the model.) The avoided 
crossings in the red circled regions reflect the (clearly small) 
effects of the orbital mixing terms on the Fermi surface struc- 
ture. The right panel shows the intercomponent Josephson 
coupling, J7b, defined in Eq. [3] as a function of the interor- 
bital coupling energy, St, computed by numerically perform- 
ing the integral in Eqs. [14] and IA7I in Appendix IA II with 
A /t = 0.002,0.004,0.008 and 77 = 0; the result confirms an 
approximate linear dependence of Jo on |5t|A , as in Eq |17l 



C. Spin-orbit term, To 

The xy orbital does not play a direct role in our analy- 
sis, but does contribute to the intra-orbital d-vector lock- 
ing r. While we will not consider the superconductivity 
in this orbital, it does affect the pairing interaction in the 
quasi-lD orbitals. It is due to the spin-orbit coupling in- 
volving the xy orbitals that we have anisotropy in the 
spin channel of the pairing interaction. More specifically, 
in both quasi- ID orbitals, the effective pairing interac- 
tion for d || z will be stronger than the effective pairing 
interaction for d _L z. This difference in the pairing in- 
teraction can be estimated as SV/V cx (r]/t) 2 as it is due 
to the electrons having intermediate states in the xy or- 
bitals (see AppdixlBlfor derivation). Now, from the BCS 
self-consistency condition log(£/A) cx t/V, the pairing 
gap will change by 



SA tSV 
-A*V^ 



loi 



t 

A 



(18) 



for small SV. Therefore, the condensation energy for the 
case d\\ z should be slightly larger by 



r = ^N(0)5(A 2 ) cx 



r/\2 t 
V l0g A 



N{0)A 2 (19) 



approximation presented in Appendix IA 11 As promised, 
Jo is paramctrically smaller than the condensation en- 
ergy per unit volume, iV(0)A§/2. 



A few points are worth noting. Firstly, this result 
is non-analytic in St, from which one concludes that it 
is non-perturbativc. In Appendix IA 21 we compute Jo 
perturbatively in powers of St, in which limit we obtain 
Jo ~ N(0)(St) 4 1 (iAo), which is an analytic function of 
St, but non-analytic in Aq. This reflects the fact that 
the perturbative expression is valid only for \St\ <C Ao, 
a physically unreasonable restriction. The perturbative 
expression does, however, match smoothly to the non- 
perturbative one when \St\ ~ A . The origin of Eq. [17] 
can be understood intuitively as arising from the quasi 
ID character of the bands. The contribution to the con- 
densation energy from the portions of the bands away 
from the crossing points (enclosed by red circles in Fig. 
I17[) is largely insensitive to orbital mixing - only in a 
neighborhood of width ~ \St/t\ about the crossing points 
is orbital mixing significant, but there it makes an 0(1) 
change in the condensation energy. These considerations 
lead to the proposed scaling expression. On the other 
hand, if the chiral pairing originated from the nearly cir- 
cular 2D Fermi surface, changes in the relative phase of 
the p x and the p y components of the pair-field affect the 
pairing gap magnitude over most of the Fermi surface, 
so Jo (or more properly, £q) must be order of the total 
condensation energ y 18 i 24 i 25 . 



than that of the case d _L z; this is the origin of the intra- 
orbital d-vector locking in Eq[2] We emphasize that the 
energy scale of this locking is much smaller than that 
of the condensation energy, which is consistent with the 
c-axis Knight shift experimen t 26 . 27 However, this small 
d- vector locking, unlike the small J, is not a characteristic 
unique to the quasi-lD model, as the pairing interaction 
anisotropy due to the same physics can also occur in the 
2D orbital. 



V. NORMAL MODES 

In the usual fashion, we can obtain an understanding 
of the low energy collective modes from the equations 
of motion derived from the non-linear sigma model 2 ^ in 
Eq. [2] The d- vector dynamics are those of coupled pen- 
dulums, while the relative phase, 0_ executes familiar 
Josephson oscillations. Looking at these modes in the 
limit that k > (spatially homogeneous modes) and 
for small amplitude deviations from the A-phase ground- 



state, in which 0_ = and d x = d y 



we find 



1 Ac Q 2 1 



^ Xsp . 



Sd x 









To + 4 Jo -4 Jo 

-4 Jo r + 4 Jo 



3d 
5d 



(20) 



where i5d a ■ e z — 0. The above equation of motion tells 
is that the small value of J in the quasi-lD model gives 
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us more soft collective modes than the 2D model. From 
this we can deduce the the gaps (or "masses") of three 
distinct normal modes, all of which vanish in the limit 
5H — > 0, i.e. at the multcritical point: 

m c = S^/Jq/xc = (7c)A 0j 
m s+ = ^2T /xsp = (7s+)A , (21) 
m s - = ^2{r + 8J )/xs P = (7.s-)A , , 

where for the relative phase mode, 7 C ~ |<ft|/i, for the in- 
phase spin-wave mode, 7 S+ ~ (ij/t) 2 , while for the rela- 
tive d- vector orientation mode, 7 S _ ~ the larger of 7 C and 
7s_. We also see, as is natural, that the in-phase spin- 
wave mode is unaffected by the interorbital coupling. 

VI. DISCUSSION 

Most studies to date have worked with the assumption 
that 7 is the active band. From an experimental per- 
spective, heat capacity measurements^ showed that the 
fraction of density of states at the Fermi level that was de- 
pleted at T c is consistent with the contribution from the 
7 band. However, the balance is delicate, as the contribu- 
tion from {a, /?} to the total density of states is similar to 
the contribution from 7. From a theoretical perspective, 
asymptotically exact calculations in the weak-coupling 
limit involving all three Fermi surfaces lead inevitably to 
the conclusion that the dominant pairing strength occurs 
among the ID bands^ 3 -. However, when stronger electron 
correlations are present, the validity of these results is 
unclear. In this regard, it is important to consider ex- 
perimental signatures that may help to discriminate be- 
tween the two possibilities for the active orbitals. This 
has been the primary motivation for carrying out the 
present analysis. 

We have shown that if the chiral p-wave superconduc- 
tivity in Sr2RuC>4 arises from the quasi-lD bands, this 
implies that the collective properties are controlled by the 
existence of a nearby multicritical point at which there 
is an enlarge emergent order parameter symmetry and 
correspondingly a set of anomalously soft "almost Gold- 
stone" soft collective modes. 

Because of the relatively weak spin-orbit coupling in 
Si'2Ru04, the in-phasc spin- wave mode is expected to 
have energy small compared to the superconducting gap. 
An estimate of the in-phase spin-wave gap can be ob- 
tained on the basis of experiments on the the c-axis 
Knight shfit=2., from which it follows that 7 S + = m s+ / Ao 
is less than 0.01. This result implies an extremely small 
value of r in Eq. [5J but does not distinguish be- 
tween different microscopic origins of the pairing. We 
have therefore focused, in particular, on the gap (mass) 



of the relative phase and relative d-vector oricntational 
modes, which are analogues of the "clapping modes" fa- 
miliar from studies of 3 He- A 24 > 25 and also investigated in 
the chiral d-wave superconductor—. The corresponding 
mode frequencies have been computed in the context of 
the quasi-2D bands, leading to the prediction that they 
would have an energy V2A (j c = 7 S _ = \/2) 18 ' 19 - We 
recover their result from our NLSM if we extrapolate the 
results to the case of strong inter-component Joscphson 
coupling, where j7o is comparable to the condensation 
energy density. However, because inter-orbital mixing is 
relatively weak, the corresponding modes are expected 
to have parametrically lower energy if the superconduc- 
tivity arises in the quasi ID bands - hence Eq. [2TJ (For 
more on this correspondence see Appendix [Pi) 

The various collective modes in Sr2Ru04 can, in princi- 
ple be detected using methods that have already been dis- 
cussed in the literature for the 2D model, which includes, 
among others, electron spin resonanc o 30 ' 31 (see Appendix 
ICj) . ultrasound attenuation and Raman scattering^. 
However, we can expect the inter-orbital 'nearly Gold- 
stone' modes to have much lower energy if the pairing 
originates primarily on the quasi ID bands, than if it 
arises on the quasi 2D band. 

Our analysis has potential implications for various ex- 
isting experimental puzzles concerning the properties of 
Sr2Ru04 near T c . If the phase transition exhibits mean- 
field behavior, an inescapable consequence is that the 
transition must be split by a field applied in the basal 
plane. The fact that this does not occur in Sr2Ru04 
suggests that fluctuations may play a significant role 
in the transition. The higher emergent symmetry of 
the proximate multicritical point would be an obvious 
source of anomalously strong fluctuations. In this regard 
it is worth noting that in two dimensions, a non-zero 
(Kostcrliz-Thoulcss) transition is still possible, despite 
the existence of gapless spin-wave like modes, but for 
an order parameter with a larger continuous symmetry, 
such fluctuations necessarily reduce the transition tem- 
perature to T = 0, or to a low temperature at which 
explicit symmetry breaking terms, or three dimensional 
couplings cut-off these fluctuations. These issues will be 
investigated in depth in a future publication. 
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Appendix A: Derivation of inter-orbital Josephson coupling 
1. BdG formalism 

Since the sum of the quasi-particle eigenenergy gives us the total ground state energy, in obtaining the ground state 
energy a given relative phase <p between two orbitals, the most straightforward method is through calculating the 
BdG eigenenergies when the relative phase is <f>. We can compute these eigenenergies from diagonalizing the spinless 
BdG Hamiltonian, 



H 



BdG \ 





A k 


Ao sin k x 





A k 




e l4 


- Ao sin 


Ao sin k x 





— £xk 


-A k 





e-^-Ao sin k y 


-A k 





(Al) 



where our basis is (u^, u y ^, v x ^, v y ^) T . This gives us the eigenenergies of 



E- 



k±l 



1 



A k ±y/B k (4>-), 



(A2) 
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where 



=Ck + ^ k + A 2 (sin 2 k x + sin 2 k y ) + 2\l, 
B k (0_) = [(£ 2 k - 6y k ) + Ag(sin 2 k x - sin 2 fc„) 2 ] 2 + 4A 2 [(& k + £ yk ) 2 + A 2 (sin 2 k x + sin 2 k y )] 
— 8A 2 . Ag sin k x sin k y cos 0_ . 



(A3) 



The ground state energy has minima at 0_ = ±ir/2 and therefore can be approximated by Eq ll3l To see this note 
that the total ground state energy 



£ (4>-) = const - \ Y^Ek+ty-) + £ k -(<M- 



(A4) 



when differentiated by </> gives us 



A^AgSin k x smk y 



E^ 



)E k - 



.)[£ k +W>o) + £k- (4>-)Y 



(A5) 



for which we used 



-) = ±- 



2A 2 . A 2 sin k. x sin k y 



E^.)[El + {^)-El_{4,.)] 
This derivative vanishes at (f> — ±7r/2, as E^±{(p — ±7r/2)'s are even in fc^ and fcy. From the second derivative 



(A6) 



£ '(^_ = 7r/2)=2£ 



A k Ag sin 2 k x sin 2 fc y 

(£ k+ i; k _) 2 (£ k+ + £ k _) 



i 



1 



£k+£k- (£ k+ +£ k _) 5 



(A7) 



which is clearly positive we see that these extrema at = ±7r/2 are minima. Lastly, we note that £o((f>—) is 7r- 
periodic in </>_, as Eq lA3l shows that each eigenenergy is invariant under the combination of — > ir — <fi- and ir/2 
rotation in the A>space. 

We now need to evaluate Eq. (IA7I) . Since the main contributions will come from the four crossing points of the ID 
orbitals, we take the following expansion around (k x ,k y ) = (27r/3, 2-7r/3): 



\ 



v(k x + k y - 47r/3) / v 2 (k x - k y ) 



± 



A 2 



A\/ q x ± + 1 + 



A 2 
A 2 "' 



(A8) 



where u = V%t/2, A = 3|<5t|/2, A = A | sin27r/3|, and g x = v(k x + k y — 4tt/3)/2X, q y = v(k x ~k y ) /2\. We also note the 
following two points: i) min(i?±) = A and ii) when E + (E-) is at its minimum, ~ A 3> A so min(_E + +i?_) ~ A. 

This leads to the following approximation: 



^ = k'(^- = V2)4E 



A k Ag sin 2 k x sin 2 k y 



_ A 4 A 4 2A 2 

2 ^ (£ k +£ k -) 3 (£ k+ + £*-) ~ 8^ 2 ~ 



d 2 g 



(£ k +£ k _) 3 (£ k+ +£; k _) 



( 2A 4 
ir 2 \v 2 

A 4 

! 8tt 2 Aw 2 



d 2 g 



/ \ 2 A 2 



A 2 



/ rr, \ 2 A 2 



A 2 



A 4 2A 2 7T 



{q 2 + A 2 /A 2 )l((j 2 + l) 2 8tt 2 Aw 2 A 2 2 



J_^A| 
16?r t t 



(A9) 



We obtained the same dependence on parameters as in FigJ^l though the coefficient came out about an order of 
magnitude larger. 
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2. Perturbation method 



We show here that once we ignore the spin-orbit coupling, Due to the C4 symmetry, this inter-orbital Josephson 
coupling is zero for the lowest order. To see this, we note that from the second order perturbation theory 



7(1), 



E 

ks 



( A k4k s c yk,A_ k c^_ ks c^_ ks + h.c.) 



E 



yk 



= -(^) 2 E 



sin 2 k x sin 2 k y A* k;ss A yk;s 



ks 



E r 



E 



E ,■ 



E 



c.c. 



yk 



where E ak = */£ 2 k + l^akl 2 , and we used 



2£„ 



(A10) 



(All) 



this result is basically due to A x r y -s k being odd in sink x r y y 

Therefore, it is from the second order inter-orbital Josephson coupling that gives rise to the dependence of the 
energy on the relative phase and the spin state d xz>yz . Given that we have completely decoupled opposite spins, 
we only need to consider the process that tunnels two spin up-up pairs and two spin down-down pairs. This can be 
calculated from the fourth order perturbation theory: 



£ j (<£-) = E ( X kclk s c yks \ kl cl k , s c ykls X- k cl _ ks c y _ ks 



Cii . — ks A— k' .sp 



x,-k',s^V 



C y ,-k' 



h.c. 



k'k;. 



(— E xk — E y k) (~E xk — E yk — E xk > — E yk i) 



1 



-E x 



^ A k A_ k A k /A_ k ' 
1 



k'k 



A\ 

2£ , xk 



A 



yk;ss 



2E. 



yk 



a?k';ss 

2E xV 



' E yk 

^yk';ss 



2£. 



yk' 



-E X k' — Eyk> 

(l-5k'k)(l-^k',-k) 

1 



E 

ks 



(— S xk — Eyk) {—E x \t — E y k — E x w — E y k<) \—E x k — E y k —Exk' — E y k< 
[(A* k , ;ss A, k , ;ss ) 2 +c.c.] 



{E xk + E yk ) s 



which gives us 



(E X kE y k) 2 

A k /4 A^sinfc^sinfcy) 2 (<5t) 4 



(E X k + Eyk) 3 (ExkEyk) 2 



t 2 A 



(A12) 



(A13) 



Note that this result is consistent with our BdG calculation, for it is qualitatively the same as taking the first term 
of Eq. (|A7[) which is much larger than the second term in the |A| 3> \6t\ limit. 



Appendix B: d- vector locking 

To see why we need the spin involving the xy orbital to lock the ci-vector along the c-axis, we need to examine the 
spin-orbit coupling part - with the xy orbital included - of the orbital hybridization of Eq. ((BJ) : 

SHkin = ^E E ^ ab ' 0- f«' C oka C -*s'- ( B1 ) 
a,b k;ss' 

This means that when we only include the spin-orbit coupling between the ID orbitals, S z will remain a good quantum 
number while S XiV will not. Therefore, it is energetically more favorable to have an equal-spin pairing with the spin 
quantization axis along the z direction, which equivalent to the (i-vector lying in the xy plane. For the same reason, 
the spin orbit coupling between the xz and xy orbital will favor the d-vector lying in the yz plane and that between 
the yz and xy orbital the rf-vector lying in the xz plane. We therefore conclude that the d-vector is locked to the c-axis 
because the spin-orbit couplings involving the xy orbital have larger effect than the spin-orbit coupling involving only 
the ID orbitals. 

We will now show that the pairing interaction anisotropy in the normal state is proportional to r] 2 . To do so, we will 
calculate how the normal state pair-field susceptibility, which is proportional to the normal state pairing interaction. 
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To account for the effect of the spin-orbit coupling, what we will calculate is the inter-orbital susceptibility of the 
intra-orbital pairs. Since the form of the spin-orbit coupling is the same for any pair of orbitals, we can expect the 
dependence on rj to be the same. Therefore, we will only look at the triplet pair susceptibility involving one pair in 
the xz orbital and another pair on the yz orbital: 

xt S C;x-y(K *«) = J2 [ dTemT ( T r4,k, a W4 ) -k^(^)^,-k,7(0)c !/ ,k,A(0))(ia 2 d • *) a p{-i& ■ <T<7 2 ) 7 A, (B2) 

CK/37A * 







ignoring the xy orbital. In this case S z is a good quantum number, so we obtain 

Xtsc-.x-yKkJty E G ya . xa (-k,-iu n -iQ;)Gy if . tXi f(k,ibj), 



(B3) 



iu n cr=t,4- 



for the d-vector in and out of plane, respectively. When we take into account that the Green function is diagonal not 
in the orbital basis but in the band basis, we can write 

|| /i, ffl- 1 V r V fa, -k, a\fi, -k, a) (fi, -k, a\x, -k, a) (y, k, a\v, k, a) (v, k, a\x, k, a) 

Xtsc; X -yK K > ni ) -p 2^ 2^ 2^ 77. 7F> TT17 77. T777 ' 



iu n - id- k) 



iu n - 6,(k) 



g n . n ^_l^ \ - \ - (y, -k,cr|/i, -k,a)(fi, -k,a\x,-k,a) (y, k, a\v, k, o)(v, k, a\x, k, a) 

Pr^.j,t^ -«u n - ill - ?M _k ) zw„-£„(k) 



(B4) 



iu n a— ^,4- /xy 



where 1/ are band indices and £k's are the normal state eigenenergies. Using the fact that ^^(j/k|yuk) (juk|a;, k) = 
and, in an appropriate basis, (aa\pa) = (aa\pa)* , a being the orbital label, we obtain the anisotropy 

AxtSC;x-y(K ifi) =XtSC,x- y (K ity ~ x\ S C;x-y ( k ' 

= - [I (yk t |ak t) (ak t |zk f) | 2 - Re((yk | |ok t) <«k t |xk t)) 2 ] 



E 



-iuj n — ifl — £ Q (— k) -!w„ — id — £p(— k) 



«w„ - £ Q (k) zw„ - £g (k) 



(B5) 



It is the transformation matrix between the orbital and the band basis that gives rise to the rj 2 dependence. 
To see this, note that the first quantized form of the normal state Hamiltonian can be written as 



hkin + $hki 



Cx + £,y £ a - £/3 



cos p 1a sin p 
-iasmp — cos p 



(B6) 



when ignoring the spin conserving orbital hybridization, with tanp = 277/ (£ x — £ y ). This gives us (y f |a f) 
e -i7r/4 s i n (p/2) and (a f |x f) = e~ ra/4 cos(p/2), so we obtain 



|(y t l« t}<<* t Is t>| 2 - Re«y t l« t)(« t N t» 2 = 2 cos 2 (p/2) sin 2 (p/2) = 
and the pair-field anisotropy at the zero temperature and frequency is 



v */2 



(e,-e y ) 2 /4+»7 2 



(B7) 



AXtSC;x-y{k,itt = 0)| T=0 = 



17 



1 

'2[77 2 -^(k)^(k)]^(k)+^(k)]' 
I 



(B8) 



Appendix C: Experimental detection 

In detecting the collective modes, we can either try to 
see the resonant response to some effective driving force 



or other excitations of Si^RuCU decay into the collective 
mode. 

We will first discuss the resonant response. Any AC 
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field that couples linearly to our collective modes can 
serve as a driving force. One well known case is the AC 
Zeeman field^i, which couples through the Zeeman en- 
ergy term T~L Z = -(^b/S)H • J2 a Sa (where fi B is the 
Bohr magneton). To see this, note that we can approx- 
imate, for a small H, S a ps XspZ x (<9 t d a ), giving us the 
spin equations of motion 



Sd x 




. Sd y . 


+ 



x <9 f H 



T + 4J -4 J 
-4 J T + 4J 

1 
1 



Sd x 
Sd,„ 



(CI) 



which are just those of a pair of driven coupled harmonic 
oscillators. The sharp difference between the quasi-lD 
and the 2D model is that the former has a double reso- 
nance peak whereas the latter has only a single peak just 
like that of the longitudinal NMR in 3 He-A22. The charge 
analogue for this would be a uniaxial AC strain along the 
a(&)-axis. This will give rise to a chemical potential dif- 
ference between the xz, yz orbitals, we expect the energy 
cost to be proportional to uSN. This can act as a driv- 
ing force for the relative phase mode, as 5N hx c dt4'/'2 
for a small u. The order of magnitude estimate for these 
resonance frequencies are ~ 0.1T C ~ 10 GHz. 

We note that the layered structure of Sr2Ru04 will 
make these resonances easier to detect. Note that since 
we expect the ci-vectors to point at the c-axis in the equi- 
librium, the Zeeman field will be applied along the ab 
plane. Given the long penetration length (~ 152nm) for 
this in-plane field, we do not expect the Meissner screen- 
ing to be significant. Also, because the lower critical field 
is very small, we can actually induce a nearly uniform 
magnetic field. 

The interorbital collective modes can also be detected 
is through relaxation processes. For instance, the phonon 
modes involving displacement of the next-nearest neigh- 
bor Ru atoms will modulate the interorbital coupling and 
hence can decay into the interorbital collective modes. 
Another is the NQR relaxation of the Ru atoms due to 
the relative phase mode, which gives rise to oscillating 
electric quadrupole moments on the Ru atoms as it in- 
volves oscillation of Cooper pair numbers between the Ru 
d xz and the Ru d yz orbitals. These electric quadrupole 
moments will relax the Ru nuclear quadrupole moments 
through the Ru atomic spin-orbit coupling. 



Appendix D: Relation to collective modes of the 2D 
model 



We show here how our collective modes are related to 
the collective modes studied for the 2D model. For every 



branch of collective modes in the 2D model, we can find 
its symmetry equivalent in the quasi-lD model. However, 
there are energy degeneracies in the 2D model which we 
expect would be broken in the quasi- ID model. We will 
show how this gives rise to the possibility of having soft 
collective modes in the quasi-lD model. 



By generalizing the results from 3He-A2^, one can see 
that there should be 12 branches of collective modes in 
a 2D chiral p-wave superconductor s] 18 i 30 . Six of them 
involve no fluctuation of the orbital degrees of freedom, 
while the other six involves the relative fluctuation of the 
p x and p y pairings (termed 'clapping' modes). 



The spectra of the six branches involving no orbital 
fluctuation remain essentially the same regardless of 
whether we take the quasi-lD or the 2D model. They 
consist of the overall phase modes, the overall amplitude 
(Higgs) modes, the two branches of total spin modes, and 
also the two branches of condensate polarization modes. 
Regardless of the model, only the total spin modes can 
be soft; the overall amplitude and polarization modes 
have gaps equal to the pair-breaking frequency 2Ao while 
the overall phase modes are gapped by the plasmon fre- 
quency. 



On the other hand, the spectra of the p x -p y relative 
modes arc strongly affected by how close the system 
is to being rotationally invariant. There consist of six 
branches - the relative phase modes, the relative ampli- 
tude modes, the two branches of the relative spin modes, 
and the two branches the spin relative amplitude mode 
(the relative amplitude oscillation out of phase by 7r be- 
tween the spin up- up and down-down pairs). When the 
system is rotationally invariant, the relative phase modes 
and the relative amplitude modes are related by 7r/4 rota- 
tion around the c-axis due to the combined U(l) symme- 
try of the overall phase and the orbital rotation. There- 
fore, up to the spin-orbit coupling, all six branches have 
the same gap, which is calculated to be v^A o 18119 . How- 
ever, when the rotational invariance is broken, this de- 
generacy between the relative phase and the relative am- 
plitude mode is completely broken. Since the p x and the 
p y pairings are nearly decoupled in the quasi-lD model, 
the relative phase mode is nearly Goldstone while the 
relative phase mode have a gap that is almost same as 
the overall amplitude mode. Due to the same reason, the 
relative spin modes are nearly Goldstone, while the spin 
relative amplitude modes have a gap close to 2Aq. 
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